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It has become traditional to assume that the Dirac structure of the phenomeno¬ 
logical quark confinement potential is scalar 0 scalar. We use the heavy quark 
expansion of the Coulomb gauge QCD Hamiltonian and the Flux Tube model to 
demonstrate that this is true but in the effective sense only. The demonstration 
contains some surprises: confinement is actually vector 0 vector and it is non- 
perturbative mixing between ordinary and hybrid QQ states which generates the 
scalar-like spin dependent potential at order Thus the existence of hybrids 

is crucial to establishing well-known spin splitting phenomenology. Finally, the 
resolution also indicates that the gluonic degrees of freedom in a hybrid must be 
of a collective nature. 


1 Introduction 

Although it has been postulated for more than 30 years, the phenomenon 
of quark confinement remains an enigmatic feature of QCD. Quenched lattice 
gauge theory and heavy quark phenomenology indicate that the static (m^ >> 
Aqcd) long range potential should be linear with a slope of 6 « 0.18 GeV^. 
The order l/m'^ quark-antiquark long range spin-dependent (SD) structure 
is also of interest and is the subject of this paper. In a model approach, 
QQ interactions are typically derived from a nonrelativistic reduction of a 
relativistic current-current interaction 


Vcouf = J jr(x)Ar(x - y)jr(y)- (1) 

In the heavy quark, nonrelativistic limit, this reduces to 

Vconf ^{r) + VsD + • ■ • ( 2 ) 

where the ellipsis denotes the spin independent interaction at C>(1 /to^) and 
terms of I^her order in 1/m. The spin dependent interaction may be param¬ 
eterized as 
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Here e = e(r) is the static potential, r = jr, — r^j is the QQ separation and 
the Vi = Vi {r) are potentials which depend on the nonperturbative structure 
of QCD. Assummg that the Vi simply follow from Eq. (1) yields the results 
shown in Table !□. Tensor and axial vector currents have spin-dependent static 
potentials and are therefore not useful for phenomenology in this context. Even 
the derivation of these simple results is beset with ambiguities and controversy. 
The reader is referred to the excellent review by Gromes in Ref. □. Finally 
we note that covariance under Lorentz transformations leads to a constraint 
between the SD potentials: e(r) = V 2 {r) — Vi{r), known as Gromes’ relatioin. 


Table 1: Nonrelativisitic limits of Eq. (1). 


r 

er 

Vi 

V2 

1^3 

Vi 

scalar 

S 

-s 

0 

0 

0 

vector 

V 

0 

V 

V'/r - V" 

2V2p 

pseudoscalar 

0 

0 

0 

P" - P'/r 



Twenty-two years ago Schnitzefl examined splittings in P-wave heavy 
quarkonia in an attempt to distinguish the possible spin dependence of the 
confinement potential. Specihcally, he assumed a vector Goulomb short range 
potential and scalar or vector long range potentials. The ratio of P-wave split¬ 
tings is then given by 


^ M{^P2) - Mj^P,) ^ {^±A^) 

M(3Pi)-M(3Po) (^±hA)’ 

where A = 14 (1) and i? = 4 (1) for the vector (scalar) case. Experimentally, 
the splittings for the J/'i/' and T systems are 0.48 and 0.66 respectively. These 
are shown as horizontal lines in Fig. 1. The upper and lower curves are pv and 
PS respectively (these have not been convoluted with the wavefunctions). Since 
the bb wavefunction is smaller than the cc, the curve must drop with increasing 
r to fit the data. Clearly only scalar confinement meets this criterion. This 
conclusion has been supported bjtmany other methods, irysluding lattice gauge 
theorjB, Wilson loop calculationffQ, and flux tube modelw. 
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Figure 1: Scalar (s) and vector (v) spin splittings vs. the quark-antiquark radius. 


A consistent picture of a QCD-generated effective scalar confinement in¬ 
teraction appears to be emerging. It is therefore disconcerting that attempts 
to build models which are “closer” to QCD (typically relativistic models) seem 
to require vector conhnement. For example, scalar Salpeter equations do not 
have stable solution^ and it appears to be impossihl£|_tp construct a BCS-like 
vacuum of QCD when scalar confinement is assumecEEIEil. This is problematical 
if one wishes to dynamically generate constituent quark masses. Furthermore, 
attempts at modeling chiral pions will be hindered by the explicit lack of chi¬ 
ral symmetry in a scalar interaction. Finally, baryons anticonfine and colour 
singlet Tkmm-Dancoff or RPA bound states are infrared divergent in a scalar 
potentiate. 

This issue was resolved in Ref. El by first performing a Foldy-Wouthuysen 
reduction of the full Coulomb gauge Hamiltonian of QCD. This immediately 
establishes that the Dirac structure of confinement for heavy quarks is of a 
timelike-vector nature. Furthermore, the spin dependent structure of the long 
range force is determined by nonperturbative mixing of the qq states with 
hybrids. Finally, if one assumes a flux tube-like configuration of the gluonic 
component of the hybrids, a “scalar” spin-dependent potential emerges in a 
natural way. 


2 Heavy Quark Expansion of Hqcd 

Our starting point is the Coulomb gauge QCD Hamiltonian 
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Hqcd = j dxi/)'l'(x) [—ia • V +/3m] ^(x)+_ffyM+Vc —5 J dxip\x)a-A{x)'ilj{x) 

(5) 

where 

Vc = Ig^ J dxdyJ-^p-ix)V-\x,y;A)Jp\y), (6) 

y“*'(x,y;A) =(x,a|(V-I?)-i(-V2)(V-I?)-i|y,6), (7) 

and 

Hym = ljd^ [77-in(x):7n(x) + B2 (x)] . (8) 

The degrees of freedom are the transverse gluon field A = A“T“, its con¬ 
jugate momentum If = n“T“, and the quark field in the Coulomb gauge. 
The Faddeev-Popov determinant is written as = Det[V • T>], with 13“^ = 
V(5“^ — being the covariant derivative in the adjoint representation, 

and the magnetic field is given by B“ = V x A“ -|-(//“^'^A*' x A'^. The static in¬ 
teraction Vc is the nonabelian analog of the Coulomb potential. It involves the 
full QCD color charge density which has both quark and gluon components, 

p“(x) = V't(x)TXx) + r’>^A\x) ■ n=(x). (9) 

The most salient feature of the Coulomb gauge Hamiltonian is that all 
of the degrees of freedom are physical. This makes it especially useful for 
identifying the physical mechanisms which drive the spin splittings in heavy 
quarkonia. 

2.1 The Foldy-Wouthuysen Transformation 

We proceed by performing a Foldy-Wouthuysen transformation on the QCD 
Hamiltonian. This is done in complete analogy to the quantum mechanical 
case where an operator is constructed which removes the interactions between 
upper and lower components of the quark wave function order by order in the 
inverse quark mass, except that the unitary transformation is now constructed 
in Fock space. The resulting Hamiltonian is given by 

Hqcd^Hfw = J dx(mqh^x)h{x) - mqxH^)xi^)) + Hym + 

VC+H 1 +H 2 + ..., (10) 
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Hi = J dxh'<{x) (D'^ - g<7-B) h(x) - {h ^ x',mg ^ rUg), (11) 

H 2 = J dxh'^{K)ga ■[E,x'D]h{x) + {h ^ xl'mq ^ mg). (12) 

In this expression h = {1 + f3)'ip/2 and X = (1 — (3)'4^/2 denote the upper and 
lower components of the quark wave function and correspond to the annihi¬ 
lation and creation operators of the heavy quark and antiquark respectively. 
The ellipsis denotes terms which are either of 0(1 /to^) or are spin-independent 
at order 1/m^. Finally D = zV -I- gA is the covariant derivative in the funda¬ 
mental representation. 

2.2 The Static Potential 

To leading order in the quark mass the Hamiltonian describes two static, non¬ 
interacting quarks. At 0{imP) the Hamiltonian reduces to Hq = Hym + Vc- 
The eigenstates of Hq may be labeled by the quark and antiquark coordinates 
and by an index which classifies the adiabatic state of the gluonic degrees of 
freedom, 

HQ\nr-,YgYq) = en{r)\nr]VgYg). (13) 

Note that we have made explicit the dependence of the gluonic degrees of 
freedom on the position of the quarks, r. 

The corresponding eigenenergies, e„(r) may be identified with the Wil¬ 
son loop potentials calculated on the lattice. Thus for example, eo(r) is the 
Coulomb plus linear potential seen long ago. Static hybrid states are collec¬ 
tively denoted \nr]YqYq) with Ur 7 ^ 0. In recent studies llJ (see Fig. 2) the 
lowest lying adiabatic hybrid potential, ei(r) has been evaluated. 

While both Hym and Vc may contribute to the linearly rising potential 
energy seen on the lattice, it is clear that the quarks may only interact with 
the flux tube via the nonabelian Coulomb interaction. Therefore the Dirac 
structure of confinement corresponds to 70 C) 7 o from the product of color 
charge densities (see Eqs. (^ and (||)). As stressed in the Introduction, this 
appears to be at odds with 20 years of quark model phenomenology. Since the 
phenomenology is based on spin splittings, it will be instructive to examine 
the 1 /m^ perturbative corrections to the static potential. 

2.3 Spin-dependent Potentials 

The spin-dependent first order correction to the static potential is given by 
(5en^(r) = {nr;YqYq\H 2 \nr;YqYq). The order 1/m term is not considered be¬ 
cause D 2 is not spin-dependent and the matrix element of cr • B vanishes. 
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Figure 2: eo(-R) (squares) and ei{R) (circles and diamonds). From Ref. 14. 


Evaluating this yields the standard classical plus Thomas precession spin orbit 
interaction 


( 1 ) ^ / ag-Lg 

\ 4m2 


u'q ' ^1? \ ^ 
Arrig J r dr 


(14) 


Thus the first term in Eq. (^, generalized to any adiabatic potential and 
therefore true for both ordinary and hybrid QQ states, is reproduced. Further 
SD corrections occur at second order in perturbation theory in Hi. This is 
shown diagramatically in Fig. 3. As indicated in the figure, Hi contains two 
operators, and a ■ B which may act on quark or antiquark lines. The 
shaded blobs represent mixing via intermediate states which must be heavy 
quark hybrids. In Fig. 3a tr • B acts on a single quark line and therefore does 
not yield a spin-dependent contribution to Vsd- Fig- 3b, alternatively, yields 
V 3 and V 4 ; while Figs. 3c and 3d give Vi and V 2 respectively. 

Thus, the application of the Foldy-Wouthuysen transformation to the 
Coulomb gauge QCD Hamiltonian shows that the spin-dependent effective 
potential may be simply interpreted as nonperturbative mixing with hybrid 
states. This makes it clear that it is possible for nonperturbative dynami- 
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Figure 3: Contributions to the second order energy shift. 


cal physics to generate an effective spin-dependent interaction which mimics 
a scalar interaction. Actually demonstrating this requires that the matrix 
elements be evaluated. 

3 Model Evaluation of the Spin-dependent Potentials 

Before proceeding to a model evaluation of the matrix elements we note that 
it is possible to make some general statements on their expected properties. 
The results of lattice gauge theory make it clear that a flux tube-like configu¬ 
ration of glue exists between static quarks. If one thinks of this as a localized 
object with an infinite number of degrees of freedom, then it is apparent that 
V 2 must evaluate to zero. This is because the electric field operator creates 
a local excitation in the flux tube at position (see Fig. 3d). This must 
then be de-excited at by the magnetic field operator. However, the two 
operators become decorrelated because infinitely many degrees of freedom in¬ 
tervene. Similarly, the long range portions of V 3 and V 4 both vanish. Thus, 
by Gromes’ relation, the only nonzero long range interaction must be given 
by Vi = —e. This is precisely the situation required for “scalar” confinement. 
It is therefore entirely plausible that an effective scalar confinement is gener¬ 
ated by nonperturbative mixing with hybrids. Furthermore, the structure of 
the spin-dependent terms depends crucially on the nature of the ground state 
gluonic degrees of freedom and clearly favors a collective rather than a single 
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particle picture of them. 

These simple expectations are borne out in an_explicit model calculatioinJ 
where the Flux Tube Model of Isgur and Paton Ej was used to evaluat£. the 
relevant matrix elements. The interested reader is referred to Ref. E3 for 
details. Evaluating the spin dependent potentials in the Flux Tube Model 
requires explicit expressions for the electric and magnetic fields in a flux tube. 
It was therefore necessary to extend the Flux Tube_Model somewhat. Here I 
take the opportunity to correct some typos in Ref. E3. 

The commutation relation between the electric and magnetic fields is 


[F;“(x), Rj(y)] = + Oig), (15) 

which implies that the electric and magnetic field operators may be defined in 
terms of the transverse string coordinate x 

~ 4(^a(« + 1)-XaH) (16) 

and 


Bl{n) 


—i d 
na i9x|(n) 


(17) 


Here k = a^/bo and the of the commutator in Eq. (^^ is taken into account 
by the index A, Xx(.n) = Xyin) and XyM = -Xx{n). 

This corrects several sign errors in Ref. 113. Note also that Eq. 38 of that paper 
should be in terms of A not A. None of the results are changed. 

Finally, evaluating the matrix elements from the perturbative expansion 
yields 


Vi(r) = -|^CFr (18) 

and 

2 

V 2 ir) = lim (19) 

N^oo 2a^ N 

In the strong coupling limit one has b = /{2a?')Cp so that the anticipated 

expression for Vi emerges in a natural way. Furthermore, V 2 approaches zero 
like 1/fV; this is also true for V 3 and V 4 . The latter point is illustrated in Fig. 
4 where the correlation of electric and magnetic fields versus separation along 
the flux tube is shown. As expected the fields become completely decorrelated 
as the number of intervening degrees of freedom becomes large. Notice that 




Figure 4: Field Correlation Functions (for N=100). 


this implies that a constituent gluon model of hybrids would not have been 
able to produce an effective scalar interaction. 

4 Conclusions 

Spin splittings and lattice calculations indicate that confinement is scalar in 
nature. This conflicts with many relativistic models of QCD which require 
vector confinement. For example, a chirally symmetric interaction is needed 
if pseudo-Goldstone pions and spontaneous chiral symmetry breaking are to 
be generated dynamically. Furthermore it appears to be impossible to build 
a stable vacuum with a scalar kernel. We have examined this issue with the 
heavy quark limit of the Coulomb gauge QCD Hamiltonian. This approach is 
physically intuitive and is simpler to interpret and implement than methods 
based on the Wilson loop. We found that the static confinement potential must 
indeed be a Dirac timelike vector. Effective scalar interactions are generated 
at order 1 /by nonperturbative mixing with hybrid states. 

We have argued that the long range spin-spin (V 3 and V 4 ) and the vector¬ 
like spin-orbit potentials (V 2 ) should all be zero since they involve field correla¬ 
tion functions evaluated between quark and antiquark. This statement follows 
by assuming that the gluonic degrees of freedom collapse into a flux tube-like 
configuration, as shown by lattice gauge theory. Alternatively, the scalar-like 
spin-orbit potential (Vi) is proportional to the matrix element of the electric 
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and magnetic fields evaluated at the same point and hence is expected to be 
nonzero. Explicit calculations of the relevant matrix elements were carried out 
in the Flux Tube Model. The model was extended to include color degrees of 
freedom and to map the chromoelectric and chromomagnetic fields to flux tube 
phonon operators. The results obtained were in agreement with our general 
arguments and with Gromes’ relation. 

A consistent picture of the Dirac structure of confinement has emerged. 
The static central potential is timelike vector while the spin-dependent struc¬ 
ture mimics the nonrelativistic reduction of an effective scalar interaction. This 
implies that it is incorrect to employ a scalar confinement kernel when doing 
calculations with light quarks. Note however that it would be acceptable to 
use scalar confinement when working explicitly in the chiral symmetry broken 
phase, i.e., with constituent quarks in the nonrelativistic limit. The work pre¬ 
sented here also implies that a constituent gluon picture of hybrids will yield 
incorrect results for certain observables. For example, Vi and V 2 would be of 
comparable magnitude in a constituent gluon model. In general, these types 
of models must fail when nonlocal properties of the gluonic configuration are 
considered. Alternatively, it is possible that they perform quite well when 
evaluating global properties of gluonics such as the hybrid spectrum. 
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